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Abstract 

A realistic study of the physics reach of a ti threshold scan at a 
future e~^e~ linear collider is presented. The results obtained take into 
account experimental and, to a large extent, theoretical systematic er- 
rors, as well as beam effects. Because of the large correlations between 
the physical parameters that can be extracted from the threshold scan, 
a multi-parameter fit is seen as mandatory. It is shown that the top 
mass, the top width and as{Mz) can be extracted simultaneously with 
uncertainties around 20 MeV, 30 MeV and 0.0012, respectively, while 
the top Yukawa coupling can be measured, with the previous three pa- 
rameters, to an uncertainty of about 35%, after assuming an external 
prior on of ±0.001. 

1 Introduction 

The study of the top threshold scan observables in a e~^e~ linear collider 
environment was pointed out long time ago |T| as a potentially high precision 
strategy for the determination of the top mass and eventually other relevant 
parameters such as the strong coupling constant, the top quark width and 
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the top Yukawa coupling. Since then many studies with increasing levels 
of complexity have been performed to obtain a quantitative estimate of the 
attainable accuracy. 

Earlier linear collider studies of the top threshold focused on the deter- 
mination of the top quark mass ^, ^, |^. A strong correlation between 
the top mass and the strong coupling constant, as{Mz), was noticed, so that 
both quantities had to be measured at the same time, through a simultaneous 
two-parameter fit @, The correlation limited the experimental precision 
that could be achieved for the top mass to about 300 MeV. On top of it, when 
next-to-next-to-leading corrections to the it cross section at threshold where 
computed 0] , they were found to be large and to perturb the determination 
of the mass at the level of about 500 MeV. 

In 1999 there was a substantial breakthrough when two new definitions of 
top mass ("potential subtracted" and "IS" P] ) were proved to be much 
less sensitive to higher order corrections than the pole mass used previously. 
As a welcome side effect, correlations between and these new masses were 
found to be much reduced , so that a determination of rrit with less than 



100 MeV experimental error and about 100 — 150 MeV theoretical uncertainty 
became feasible The decrease in the correlation can be understood from 
the fact that the most sensitive information comes from the position of the 
IS* resonance peak (although severely smeared after initial state radiation 
(ISR) and beam effects), which is at an energy £'15 = 2mt — V^t(as), where 
rrit stands for the top pole mass, used previously, and V^f(as) stands for the 
tt binding potential, and depends almost linearly on a^, therefore heavily 
correlating both parameters. 

In order to help disentangle these two variables, the top momentum dis- 
tribution due to the top Fermi motion was advocated already several years 



ago [|Tj, |T2[. Simulation studies showed that the peak position of that distri- 



bution was quite robust against ISR and beam effects and changed linearly 
with the top mass while being insensitive to the strong coupling constant, 
providing, therefore, an additional handle for the disentangling of both vari- 
ables 0, |. 

Finally, in addition, the use of the top forward-backward charge asym- 
metry was suggested in the past as a way of getting direct information on 
the top quark width p^ , p!3[] . The large width of the top quark is responsi- 



ble for the overlap between the IS" and IP states whose interference causes 
a forward-backward charge asymmetry. Simulation studies showed that the 
measurement was feasible, although the attainable accuracy was rather lim- 
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ited i. 

In all these observables some modest sensitivity to the influence of the 
top Yukawa coupling, entering the ti potential was expected. Again, Monte 
Carlo studies showed that indeed the sensitivity was quite low 

For the top width and the top Yukawa coupling studies, the approach 
followed so far was a single parameter determination assuming no relevant 
uncertainty in the other parameters (top mass and strong coupling constant). 

The study presented in this paper continues the work of ref . flOl , extend- 
ing it and completing it. It completes it because it includes not just the 
cross section but also the other two relevant observables, momentum distri- 
bution and forward-backward asymmetry. It extends it because it explores 
the feasibility of extracting information about all the relevant input param- 
eters simultaneously, that is, not just the top mass and the strong coupling 
constant but also the top width and the top Yukawa coupling. 

The outline of the paper is as follows: in section 2 the conditions in which 
the present analysis has been developed are summarized. In section 3 the 
results for a "standard" two parameter fit to the top mass and the strong 
coupling constant are presented and the sources of the correlations obtained 
are discussed. The case for a multi-parameter approach is then presented in 
section 4 and the technical solution is described. Using this multi-parameter 
approach, section 5 deals with the discussion on the actual sensitivity to the 
top quark width and section 6 with the prospects for a measurement of the 
top Yukawa coupling. Finally, section 7 summarizes the conclusions of this 
study. 



2 Simulation Input 



For the present study, the simulation conditions have been assumed to be 
identical to the ones used in namely, the TESLA beam conditions of 



ref. [|T5| were assumed and the detector effects were simulated using the 
SIMDET algorithm as described in ref. 
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The tt observables (cross section, top momentum distribution and forward- 



backward charge asymmetry) have been computed using the TOPPIK code JTT 
T3| , 1^, IS] including the latest theoretical predictions as discussed in ref. [IS]^ 



very recent update Q of ref. ||l9[ results in shifts in the cross section of up to 1.5%, 
weU within the assumed theoretical error. The new spin-independent \/m? potentials 
in pol can be seen to agree with those obtained by Pineda in while those in 119] 
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For the studies presented in this work, the IS* mass definition has been 
used. The actual values of the input parameters used in the calculation 
are mt{lS) = 175 GeV, a,{Mz) = 0.120, Mh = 120 GeV and the top width 
and Yukawa coupling as predicted in the Minimal Standard Model. For the 
experimental selection studies, the signal and the relevant backgrounds have 



been generated using PYTHIA p2|. The event simulation is discussed in detail 



in ref. \23 



For the cross section analysis, purely hadronic decays of the tt system, 
together with events in which only one of the top quarks decays hadronically 
are used. This results in an event selection efficiency of 41.2 % (over the 
complete ti sample) with an estimated systematic uncertainty in the selection 
efficiency of 3% and a remaining background cross section of about 0.0085 
pb. The same sample is used for the study of the top momentum distribution 
and, hence, the efficiency remains the same, with an estimated systematic 
uncertainty in the peak value of the momentum distribution of about 4%. 
For the forward backward asymmetry measurement, only the sample in which 
one of the top decays semileptonically to an electron or a muon, allowing to 
tag easily the top charge, is used. This results in an event selection efficiency 
of about 14 % with negligible background and systematic uncertainties. 

For the first time, theoretical uncertainties have been included in the fit. 
Following ref. []TU| a 3% uncertainty in the total cross section, common to all 



center-of-mass energy points, has been assumed. This has to be considered 
as only a first attempt at quantifying the infiuence of the theoretical error 
in the results No estimate of theoretical systematics is available for either 
the top momentum distribution or the forward-backward charge asymmetry. 
However, as it will be shown below, these two observable have a rather limited 
weight in the final results. 

The threshold scan has been assumed to consist of a luminosity of 300 fb~^ 
uniformly distributed in 10 scan points: one of them well below threshold for 
a direct background determination and the other 9 distributed symmetrically 
around the mt(lS') value with a 1 GeV spacing between each other i. Since 



disagreed. 

^For instance, we have noticed that assuming, instead, that the 3% theoretical error is 
totally uncorrelated between energy points (which seems rather unlikely), the effect of the 
theoretical error becomes much more prominent. 

^Some studies carried out in the past using just the cross section and fitting only the 
top mass and the strong coupling constant showed already that sonic modifications of 
the scanning strategy could allow to decrease some parameter errors [lO|. However, when 
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the top mass will only be known with a moderate precision at the time 
the top threshold scan starts at a linear collider, it will be impossible to 
choose the center-of-mass energy points to lie precisely at the values we have 
chosen relative to 2mt. However, it has been checked that assuming a prior 
precision on of around 500 MeV (coming from measurements at LHC or 
in the continuum region at a linear collider), the effect of not choosing the 
optimal values for ^/s is very small. 

The expectations for the three observables studied here are shown in 
fig. in together with the corresponding expected experimental errors. It is 
clear from that figure that, while the measurement of cross section will be 
very precise, the peak of the momentum distribution will be determined with 
a moderate precision and the forward-backward asymmetry with a rather low 
precision. These uncertainties should be kept in mind when analyzing the 
sensitivity of each observable to the input parameters in the next sections. 
Both things put together will provide a feeling for which measurements are 
actually important for the determination of each one of the parameters dis- 
cussed. 



3 The Top Mass and the Strong Couphng 
Constant 

To start with, a two parameter fit, with rrit and as{Mz), is performed, as 



in ref. [110| but now using the larger integrated luminosity mentioned in the 
previous section and including also the information in the forward-backward 
charge asymmetry (Afb) and position of the peak of the top momentum 
distribution, as outlined also in section 0. 

The sensitivity of each one of the observables to the top mass can be 
gleaned from fig. ^ while the sensitivity to the strong coupling constant can 
be seen in fig. |[ Given the scale of the experimental uncertainties shown in 
the previous section, it is easy to see that no relevant information in these 
two parameters can be expected from the forward-backward asymmetry. The 
peak of the momentum distribution is fairly sensitive to the top mass while 
rather insensitive to and, even within the modest experimental precision, 
should provide valuable information on m^. Finally, the cross section is very 



dealing with four parameter fits, as we are doing here, the optimization of the scan strategy 
is not so obvious and for the moment it has not been attempted. 
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Figure 1: The expected scan results for the cross section, the peak of the top 
momentum distribution and the forward-backward charge asymmetry in the 
conditions described in the text. 
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sensitive to both the top mass and as but, since we are using the IS' top mass 
definition, while most of the sensitivity to the top mass is in the (smeared) 
threshold position, the sensitivity to comes mainly from the cross section 
above threshold. This indicates already that the top mass determination 
would benefit from investing the luminosity of the scan points which are above 
threshold in the threshold rising slope instead. The resulting uncertainties, 
including only experimental errors, are the following: 

Amt = 16 MeV Aa, = 0.0011 p = 0.36 , (1) 

where p is the correlation coefficient between rrit and a.^- If the 3% theoretical 
normalization error in the cross section is included, the results change to: 

Amt = 16 MeV Aa, = 0.0012 p = 0.33 , (2) 

As it can be seen, the change is extremely small. In the rest of the paper, 
unless explicitly said otherwise, all numbers and figures given will include the 
effect of the theoretical error. Figure ^ shows the correlation plot between 
the top mass and a<j resulting from the two-parameter fit- 

To quantify the actual contribution of the different observables to this 
two-parameter fit, it has been repeated using only the cross section. In this 
case the results are 

Amt = 25 MeV Aa, = 0.0019 p = 0.76 , (3) 

The difference between this result and the one above has been traced back 
to the contribution of the peak of the momentum distribution, as expected. 
The sensitivity of the momentum distribution is such that it gives raise to 
a negative correlation between the extracted values of mt and which, 
when combined with the positive correlation coming from the cross section 
measurement, leads to substantial reductions in the overall errors of both mt 
and as- 

Once the possibility of measuring precisely the top mass and as with a it 
threshold scan has been established one may turn the attention to measuring 
other quantities, like the top quark width and the top quark Yukawa coupling. 

4 Multi-parameter fit strategy 

So far, only single-parameter (or at most two-parameter) fits have been tried 
to the threshold scan observables to study the potential for the determination 
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Figure 2: Sensitivity of the observables to the top mass. The different mark- 
ers correspond to Amj = 200 MeV intervals. 
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Figure 3: Sensitivity of the observables to the strong couphng constant. The 
different markers correspond to Aag = 0.004 intervals. 
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Figure 4: Ax^ = 1 contour as a function of mt{lS) and as{Mz). 
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of the top width and Yukawa couphng. Nevertheless, as we have seen, the 
expected experimental error in each one of the observables strongly suggests 
that very likely the cross section measurement will dominate the parameter 
determination (with only a small improvement eventually expectable from 
adding the information from the forward-backward asymmetry and the top 
momentum distributions). 

Therefore, as already stressed above, non- negligible correlations might be 
expected between the four parameters (top mass, strong coupling constant, 
top width and top Yukawa coupling) and, in this scenario, the only valid 
approach is a multi-parameter fit strategy. 

For that, we need predictions of the three threshold observables studied 
in this work as a function of the four free parameters. In spite of the fact that 
the TOPPIK code speed has sizably improved over the previous versions, the 
time needed to run it after the convolution of the predictions with ISR and 
with the beam spectrum makes its use for a 4P fit impossible. To cure this 
problem and produce the necessary predictions in an affordable amount of 
time we have used a multidimensional interpolating routine based upon the 



algorithms of ref. |]2^. We have checked that, within the parameter intervals 
relevant for the fits discussed in this work, that interpolation produces results 
which reproduce the exact predictions with the required accuracy. 



5 The Top Quark Width 

Earlier attempts have been made in studying the determination of the top 
quark width (Ft) from the it threshold scan [^, 
less than perfect agreement with each other, at least apparently. 

Figure |^ shows the sensitivity of the threshold observables to the top 
width. As it can be seen from the figure, there is sizable sensitivity both at the 
peak of the cross section and at lower center-of-mass energies. For low values 
of the top width, the peak structure of the IS* resonance becomes apparent 
while for large values it disappears. The other two observables {ApB and peak 
of momentum distribution) have an even larger sensitivity, which is enhanced 
in the energy points above threshold. Unfortunately, as already stressed, the 
accuracy of their experimental determination is much poorer. Because of 
that, in the scenario presented in this work, the determination of the top 
width can be expected to be dominated by the cross section measurement. 

A three-parameter fit, with rrit, as and F^ leads to the following overall 



14, 25, 



with results in 
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Figure 5: Sensitivity of the observables to the top width. The different 
markers correspond to AFj = 400 MeV intervals. 
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uncertainties: 



Amt = 19 MeV Aa, = 0.0012 ATt = 32 MeV , (4) 

with all correlations between the three parameters being below 50%. The 
32 MeV uncertainty on the top width corresponds to about a 2% measure- 
ment. This has to be compared with the 18% uncertainty reported on ref. [^]. 
Several factors account for the decrease in the reported error: 

• The integrated luminosity assumed now is six times larger than that 
assumed in (30 fb~^ instead of 5 fb~^ per point). This accounts for 
a factor ~ 2.4. 

• The selection efficiency for top-antitop events has been taken from the 
work in |p3i, and it is about 41% while the one used in [HI was of about 



25%. This accounts for a factor ~ 1.3. 

• The present studies are based on a newer version of the TESLA ma- 
chine. As a consequence the beam spectrum is now sharper than pre- 
viously assumed. This results in an improvement of the top width 
determination by a factor ~ 1.5. 

• A substantial part of the sensitivity in the center-of-mass energies below 
the maximum of the cross section was lost in the scan of [Q, because, 
although the scan was centered around 2mt, the position of the peak 
in the cross section was shifted by about 2 GeV toward lower y/s, 
because of the pole mass definition used for mj. In contrast, both PS 
and IS masses result, essentially by definition, in a cross section with 
a maximum at 2mt, so that the scanning strategy used here catches 
those energies below the maximum with substantial sensitivity to the 
top width. This difference can be clearly seen in figure |^. This accounts 
for a factor ~ 1.8. 

When all these changes are taken into account together, a factor ~ 8.5 
difference is obtained and therefore good compatibility is found between the 
results presented here and those in [Q. 

To disentangle the contribution of the different observables to the top 
width determination, the fit has been repeated using only the cross section. 
The results are 

Amt = 34 MeV Aa, = 0.0023 ATt = 42 MeV , (5) 
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with correlations as large as 80% between the top mass and as- The differ- 
ence between this fit and the result above can be traced back completely to 
the contribution of the peak of the momentum distribution in determining 
mt. The contribution of the forward-backward asymmetry, introduced be- 
cause conceptually should be the cleanest observable to see the effect of the 
top width due to the overlap between the IS" and IP states is in practice 
negligible. 

It is important to stress here that, being the top quark so heavy, a 2% 
determination of the top quark width can be very useful in constraining mod- 
els of new physics which would predict new particles that could be produced 
on top quark decays. The precise determination of the top width from the 
threshold scan allows to put constraints which are independent of the final 
state particles produced. 



6 The Top Yukawa Coupling 

Measuring the top Yukawa coupling could provide an important test of the 
Higgs mechanism for generating fermion masses. The exchange of a Higgs 
boson between the top and anti-top produced at threshold has been taken 
into account in the theory prediction by adding a Yukawa potential to the 
QCD it potential The modified potential can have measurable effects 
in the observables studied here. However, figure |^ shows that the sensitivity 
of the total cross section to the Yukawa coupling is not very large and is not 
better in the forward-backward asymmetry while is unexisting in the peak 
of the momentum distribution. Following the same arguments given in the 
previous sections, we can expect that the Yukawa coupling determination 
will be completely dominated by the cross section measurement as well. 

As an unrealistic starting point, a one-parameter fit is performed, fixing 
all parameters except for the top Yukawa coupling, A^. The fit returns an 
asymmetric uncertainty: 

_+0.18 (p\ 
— -0.25 • l^OJ 

Given the lack of sensitivity, one can try to see whether there could be any 
gain obtained by relaxing somewhat the assumptions concerning systematic 
errors. In particular, the systematic error in the cross section determina- 
tion has been lowered from 3% (taken from |2^) to 1%, which seems like 
a reasonable educated guess, given the level of understanding achieved at 
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electron-positron machines like LEP, where selection systematics routinely 
achieved the few per mil level. Also, maybe with less justification, the the- 
oretical error in the overall normalization has been lowered from 3% to 1%. 
Assuming the 1% errors, the uncertainty in the one-parameter fit decreases 
to: 



]t _+0.14 



/^From now on, these lower systematic errors in selection and normalization 
will be assumed on all fits. 

The next step consists on leaving the top mass and as free in the fit while 
fixing the top width to its Standard Model value and including an external 
constraint on as{Mz) of ±0.001. The constraint should come from a different 
determination of as, like the one available, for instance, at GigaZ ||2^. Under 



these conditions, a three-parameter fit leads to the following precisions: 

Amt = 27 MeV Aa, = 0.001 (constraint) ^ =+[J-^^ , (8) 

A* 

with correlations large, up to 80%, particularly among mj and A^. 

Finally, one could also try to leave the top width free in the fit, and 
perform a four-parameter fit with an external constraint on as{Mz). The 
results are: 

Amt = 31 MeV Aas = 0.001 (constraint) 
AF, = 34 MeV ^ =t'oil ■ (9) 

The simultaneous determination of the four parameters is possible without 
a large increase in the resulting uncertainties. Correlations remain similar, 
with a maximum of 83%, again among rrit and A*. 

As can be seen, a realistic determination of the top Yukawa coupling, 
which has to be done simultaneously with that of the top mass, is very 
challenging, although not impossible. 



7 Summary 

For the first time a simultaneous 4-parameter (top mass, strong coupling 
constant, top width and top Yukawa coupling) fit to the expectations for 
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three tt threshold scan observables (cross section, top momentum distribution 
and top forward-backward charge asymmetry) has been carried out. 

In a complete 4-parameter fit to the expected threshold observables for 
a total 300 fb~^ luminosity scan using the TESLA machine, the parameter 
correlations obtained are very important (up to 83%), justifying the use of a 
full multidimensional approach. 

The expected experimental uncertainties in each observable, as well as 
the sensitivity of each observable to the fitting parameters, have been scru- 
tinized to understand in detail the origin of the correlations obtained in the 
parameter determination. 

The outcome of the fits shows that the prediction of a high precision 
determination of the top mass, with an experimental accuracy better that 30 
MeV and of the top width, with an accuracy at the 2% level, is quite robust. 

Measuring the top Yukawa coupling with a top threshold scan looks dif- 
ficult. Even with somewhat optimistic assumptions, the best error which 
can be expected with the above luminosity is above 30%, assuming a Higgs 
mass of 120 GeV. The situation should become significantly worse for heavier 
Higgses. 

For the first time, an estimate of the theoretical error in the cross section 
prediction (a 3% normalization error from |T^) has been included in the 
fits and has been shown to have little effect on the overall result. However, 
a more realistic theoretical error model might lead to more severe errors, 
particularly on the top mass and as{Mz)- 
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